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ABSTRACT 


Since Newton, the stumbling-block in celestial me- 
chanics has been the three-body problem. Only restricted 
cases have yielded solutions. This paper describes a new 
tool for solvingthe perturbation problem— "spinors,” a vec- 
torial concept (in the complex plane) originally used in 
atomic physics but now applied to celestial mechanics. The 
two -body problem is discussed and two devices are intro- 
duced: "pseudo-time" and the Levi-Civita transformations. 
Then spinors are shown to make possible the extension of 
these transformations to three dimensions. Perturbation 
problems are discussed; notably a body moving around the 
sun and being acted on by the orbital motion of a major 
planet. 
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ON THE APPLICATION OF SPINORS 
TO THE PROBLEMS OF CELESTIAL MECHANICS* 


by 

Karl Stumpff 

Goddard Space Flight Center 


REGULARIZATION OF THE KEPLER PROBLEM- 
TRANSFORMATIONS OF LEVI-CIVITA 

Application of spinors to the problems of celestial mechanics was first made by P. Kustaanheimo 
and E. Stiefel. Originally, spinors were the creation of atomic physicists and mathematicians; now 
astronomers have taken them up. Spinors are symbols by which spin (the rotation of an elementary 
particle) can be mathematically described. In celestial mechanics spinors are no longer connected 
to their conceptual origins; they have become a formal aid in simplifying theory and its application 
to special problems. 

Spinors may be defined in different ways. One definition identifies a spinor with an arrow in 
the complex plane leading from a complex number z x to another complex number z 2 . Mathemati- 
cally we may represent a spinor by the matrix 

s . ("' ’■) , 

\ U 2 V 2 / 

composed of the four real numbers u i? v. , where the complex numbers are z x - + ivi and 

z 2 = u 2 + iv 2 . Or we may consider the spinor S as a four-dimensional vector 


S 



♦This essay originates from lectures given at the Goddard Space Flight Center, Laboratory for Theoretical Studies, in November, 1966. 
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Mathematicians have elaborated a complete algebra of spinors (Reference 1), closely related 
to the algebra of matrices or vectors. It would take too long and lead too far from the main object 
of this paper to give even a short introduction to spinor algebra. However, we may use instead of 
spinors the corresponding matrices or vectors or pairs of complex numbers, the algebraic proper- 
ties of which are well known. 

Our object is to solve the equations of motion of celestial-bodies. This means (1) the undis- 
turbed motion (in conic sections) of planets, comets , or artificial satellites around their central 
body, and (2) their motion when perturbed by the attraction of other bodies of the solar system. The 
differential equations of motion and their solutions in classical form have unsatisfactory features: 
lack of symmetry, lack of universality, and the existence of singularities. Spinors or their vec- 
torial analogies are found very helpful in removing these difficulties. 

Consider, first, the simple case of the two-body problem in the Cartesian xy -plane. There, 
symmetrization, regularization, and generalization can be performed in an elementary manner by 
convenient transformations of the coordinates and the independent variable. (Spinors are not yet 
needed; they will be introduced later, for studying the motion of celestial bodies in space.) 

The equations for the undisturbed Kepler ian motion of a mass point m 2 around a central mass 
m 1 in the xy-plane are: 


x + ~ = 0 , y + ^ - 0 , r - Vx 2 + y 2 , (1) 

provided that the dimensions of time, mass and length are chosen so that k 2 (mi + m 2 ) - 1. The gen- 
eral solution can be simply obtained if we use (instead of time t) the eccentric anomaly E as an in- 
dependent variable. Let 


x - r cos (v - o>) , y r sin (v - a>) , 

where r and v are polar coordinates and w is an arbitrary constant. Furthermore 

r cos v - a(cosE-e), rsinv - a fl-e 2 sin E , 

where a and e are two other arbitrary constants. Finally, E may be calculated from Kepler T s 
equation 


E - e sinE = a 3/2 (t -T) . 

The time of perihelion passage T is the last of the four independent integration constants required. 
The equation of the trajectory in the r ? E plane is 


r - a( 1 - e cos E) . 
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This is the equation of an ellipse with one focus at the origin. This simple solution serves only for 
elliptic orbits where eccentricity e < 1. For hyperbolic orbits (e > 1), E is an imaginary angle and 
we must change the algorithm for calculating the coordinates as functions of the time. In the case 
of parabolic orbits (e = 1, a = ™) y the formulas lose their meaning; in fact, numerical calculations 
become intricate for orbits with eccentricities close to unity. These various difficulties have been 
mastered, in classical celestial mechanics, by using different methods of calculation for different 
types of orbits. 

The lack of universality can be overcome by using for the new independent variable the "pseudo- 
time," s = ya(E-E 0 ) , where E 0 is the value of E at an arbitrary initial time t 0 . Then 

ds dE _ X 

dt ” 1* dt " r 


is positive and real for all types of orbits, and s itself is zero for t = t 0 . In the parabolic case, 
where a ~ ^ and E = 0, 


s 


£ 


dt 

r 


remains finite; the discontinuity at e = 1 disappears if s is used as independent variable instead of 
E or t . The coordinates x, y of a moving body can be expressed as functions of s in forms that do 
not depend on the type of orbit described. 

This transformation also solves the problem of regularization. The energy-integral is given by 


where v is the velocity and a is the major semi-axis of the conic section along which the planet is 
moving. This equation shows that in the case of a collision between the two bodies, v approaches 
infinity as r approaches zero. Equations 1 show that the acceleration approaches infinity. Nu- 
merical computation of the orbital motion is difficult if the body, at its peric enter, approaches the 
central mass very closely. Introducing s instead of t as independent variable eliminates this 
singularity. If primes denote differentiation with respect to s, then, for any function q(t), 


Therefore 


/ ds 

dt 


1 , 
T q 


+ y 

w- 2 
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The velocity, expressed with respect to pseudo-time s, (given by the square root of the left side of 
Equation 3) approaches zero as r approaches zero. The formulas giving x(s), y(s), and z(s) in a 
spatial xyz -system are not only valid for all types of orbits, but they are not affected by singularities 

Pseudo-time s may be derived as a function of r = t - t 0 from the differential relation 


dr - r ds . 

The integral of this equation is analogous to Kepler’s equation; it is likewise valid for all types of 
orbits. 

This is more apparent if the equations of two-body motion are written with s as an independent 
variable. Then, from x = x'/r, second differentiation gives: 

_ x 1 /x" x ' r 

X = ' 7*“ = T \ T ” r 2 


i.e. 


_ x ' r ' - x y' r' - y 

x - * y" " r 

These equations must be completed by differential relations for r(s). From 

r 2 - x 2 + y 2 , 


(4) 


differentiating twice gives 


+ r r 


+ xx 


yy 


+ y (x' r' - x) + y (y' r ' - y) 


i.e., 


rr" = x ' 2 + y ' 2 ” r . 

Putting in Equation 3 for the square of the velocity gives 

r 1 9 _ 1 - r" 

r" = 1 - — or - : a 1 - — - — - const. 
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Therefore, r(s) is a solution of 


r" + a 2 r - 1 or r'" + a 2 r' - 0 . 


Another differentiation of Equations 4 gives 

r x'" + r ' x " = r'x"+x'r"-x' , 


i.e., 


rx'" = x' (r" - 1) , x" + ^ - — x' = 0 , x'" + a 2 x - 0 . 

Therefore, x, y , z, and r satisfy, as functions of s, differential equations of the same type: 

q" + a 2 q' = 0 , 

where a 2 is a real constant which is positive for ellipses, zero for parabolas, and negative for 
hyperbolas. 

Integration of these linear differential equations is very easy. If q is the general symbol for 
r, x, y, or z, then: for 


a 2 > 0: q - q Q + q x cos as + q 2 sin as , 


for 


a 2 - 0: q = q Q + q x s + q 2 s 2 , 


for 


a 2 < 0: q - q 0 + qj coshes + q 2 sinh /3s , (/ 3 2 = ~a 2 > o) . 


These expressions are still distinct for different types of orbits. They receive a universal form 
from the introduction of the Tl c -functions” (Reference 5, Chapter 5). Let k = as. Then is always 
real, and 


c Q (\ 2 ) - cos k , 




sin k 

* 


C 2 (\ 2 ) 


1 - c os 
\ 2 


K 


C 3 M 


k - sin k 

' k 3 (5) 
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are functions of K 2 and real for all types of orbits. In the parabolic case, ^ is exactly zero, and 
the c -functions become constants, c n = l/n! . Between the c -functions there exist the relations 


K 2 c 


■S (c. .") 


- c „-l s 


r 


s n ds - 


, n + 1 


( 6 ) 


From these relations it follows that the general solution of the differential equation 

q'" + a 2 q' = 0 


may be written 


q - q 0 + c i 5 q 0 ' + c 2 s 2 q 0 " , 

if q 0 , q 0 ' , q 0 " are the values of q, q', q" for s = 0. Differentiating these expressions three times 
(regarding the definition of the c functions and the relations between them) shows that q"' = -a 2 q' . 
Note that q 0 and q 0 ' may be arbitrarily given as initial conditions. Then q 0 " can be computed, in 
the case of r and x, from 


r 


0 


1 - a 2 r 0 , 



To derive the time r = t - t 0 as function of s, integrate 


dr = rds = (r Q + Cj sr 0 ' + c 2 s 2 r 0 ”) ds . 

Using the integral relations of Equation 6 for the c -functions gives 



From this "main equation" of the two-body problem, which is analogous to Kepler's equation, 
s = s(r) may be derived by iteration. 

The above solution of Kepler's problem is universal and free from singularities, but not unique. 
Levi-Civita has proposed another, simpler and possibly more suitable as an approximate solution 
when perturbations are present. Levi-Civita transforms the Cartesian coordinates x, y by 
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putting 


x - £ 2 - v 2 > y - 2£v > r - £ 2 + 7 j 2 . (7) 

The trajectory of the motion in the xy-plane (the positive x-axis being directed to the perihelion) is 
given by 


x - r cos v t 


y = r sin v , 


r(l + ecosv) - a (l - e 2 ) 


and is transformed into 


r + ex - (£ 2 +t? 2 ) + e(£ 2 -7? 2 ) - a (l - e 2 ) , 


i*e., 


<f 2 _ 

a(l - e) + a(l + e) ~ 1 ' 

This is the central equation of an ellipse, the semi-axes of which are /5( l - e) and /a( 1 + e), i.e. 
the square roots of the perihelion and aphelion distances of the true orbit (see Figure la). If, in 
the xy-plane the body moves around the sun in the direction APA , it moves in the £77 -plane only 
along one half of the ellipse. In other words, one revolution in the £ 77 -plane corresponds to two 
revolutions in the xy-plane. Indeed, if £ and 77 are expressed by polar coordinates R, 

£ - R cos x » V ~ Rsiny, 

From these equations and Equations 7, it can be shown that 

2 _ X + Vx 2 + y^ 

£ 2 


9 -x + Vx^+y 

-n = 2 — 


R - /r , 
v 

X ' 2 ‘ 
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See Figure lb. Also: 



a (1 - e ) — I*- 
(a) 



Figure 1 — Trajectory in Cartesian 
and Levi-Civita coordinates. 


R 2 ( 1 + e cos 2y) - p - a (l - e 2 ) - r( 1 + e cos v) . 

Differentiating Equations 7 with respect to time gives 

x = 2(^-7777), y = 2(^77 + 77^). 

Therefore 

V 2 - x 2 + y 2 - 4 (^ 2 + ?7 2 ) (^ 2 + 77 2 ) = 4 r (^ 2 + -77 2 ) . 

Hence, the new expression of the square of the velocity in 
the ^77-plane will be 



and, from the energy integral for Keplerian motion, 


it follows that 


1 1 
2 r 2 4ar 


This expression approaches infinity as r approaches zero but, if the pseudo-time s is used as in- 
dependent variable instead of t , 


+ V 2 = 4^ • (8) 

In Levi-Civita coordinates, the velocity of impact is 1/V2. The movement in the ^ 77 -plane is there- 
fore regular and finite for all types of orbits and all phases of the motion. In the case of a parabolic 
orbit (a = »), the velocity is constant = l/l^2. 

The equations of motion in the plane, 


X + 


X 



0 , 



0 , 
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now take the form 


V" + V (jt 


f ,2 + vi^ 




0 , 


0 - 


Equation 8 for 2 + V 2 gives: 


£" + w 2 £ - 0 , 7j" + OJ 2 77 - 0 , 6J 2 = 4^- • (9) 

These differential equations are of the same type as before. As +^ 2 = 0, etc., the general 

solution may be written 

£ = £ 0 + C t s£ 0 ' + c 2 s 2 f 0 " , etc. [c. = Cj (w 2 s 2 )] . 

But, as £ 0 " - £ 0 , therefore, with K = cus, 


£ = £0 t 1 C 2 ) + C 1 s ^o' • 

and, because i-\ 2 c 2 =c 0 , 

?' = £ 0 c o + Co' C 1 s • 

V - Vo c 0 + V C 1 s ■ 

[ c i = c i(" 2s2 )] (10) 

or, in the elliptic case (oj 2 > 0), 

sin &js 

f - f 0 cosa;s + <f 0 s ~ ’ 

sinws 

V - V 0 c °s cjs + i ? 0 s • ( 11 ) 


This solution is much simpler than the first one discussed. 
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The main equation is the same as before: 


r 



s + r n 


+ r n 


c. = c. (a 2 s 2 ) . 


But, as r 0 " - 1 - a 2 r 0 , 1 - c 3 a 2 s 2 = c t f and a 2 = 4cj 2 ; 

t - r Q Cj s + r Q ' c 2 s 2 + c 3 s 3 . 

Here the argument of the c -functions, however, is not s 2 but 

a 2 s 2 = 4 oj 2 s 2 . 

The Levi-Civita transformation is only valid for plane coordinates. The next section shows 
that it is possible (though difficult) to extend this method to the case of three dimensions. Here, 
spinors prove a helpful device. 


LEVI-CIVITA TRANSFORMATIONS IN SPACE-INTRODUCTION OF SPINORS 

The solutions of the plane two-body problem (given in Equations 10 and 11, expressed in Levi- 
Civita coordinates and with pseudo-time s as independent variable) prove to be simple, independent 
of the type of orbits, real, and finite even for collisions, and therefore easily computable even for 
small separations. The two-body problem lent itself to simplification through the use of the Levi- 
Civita transformations, as it could be restricted to two dimensions by the choosing of a convenient 
plane system of reference. The Levi-Civita transformations cannot be extended easily to three 
dimensions— a disadvantage in computing disturbed motions, as disturbing forces generally have 
components in all directions. Nevertheless, this extension can be made. 

P. Kustaanheimo and E. Stiefel (References 2 and 3) have succeeded in constructing Levi- 
Civita coordinates in three dimensions by using spinors. A spinor is defined by an arrow in the 
complex z-plane, leading from a complex number z l = u t + i Vj to another z 2 - u 2 + iv 2 . Consider 
such a spinor: 


S 


12 


(*i- z s) 



Then, 
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and 


are the conjugates of z x and z 2 . Let us define four values r, x, y, z by 

r + X = 2 ZjZj = 2 ( u i 2+v i 2 ) » 
r - x = 2 z 2 z 2 = 2 (u 2 2 + v 2 2 ) , 
y + iz a 2Zj z 2 = 2(u 1 u 2 + v 1 v 2 j + 2i (u 2 v j - u x v 2 ) , 

y - iz = 2z 2 Zj = 2( Ul u 2 + Vj v 2 ) - 2i(u 2 Vj -Uj v 2 ) , 

These equations give the relations 

f = Z 1 Z 1 + Z 2 Z 2 = U l 2 + V 1 + U 2 2 + V 2 ' 

X - Zj Zj - Z 2 Z 2 - U 1 + Vj 2 - U 2 2 - V 2 2 , 


y - z, 2 j + z 2 z, - 2 (uj u 2 + Vj v 2 ) , 

z = i (z 2 ¥j - z, z 2 ) = 2(u 2 Vj-Uj v 2 ) . (12) 

These four values are not independent, though being functions of the four elements of the spinor 
s 12 . Because 


(r + x) (r - x) - (y + iz) (y - iz) - r 2 - (x 2 +y 2 + z 2 ) - 4 (z, z, • z 2 z 2 - z, z 2 • z 2 z,) - 0. 


Therefore 


r 


Vx 2 + _ y 2 + z 2 ~ 


/TV+y, 2 ) 


is positive and equal to the absolute magnitude of the vector 
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From Equations 12, to every spinor S 12 = (z 1? z 2 ) there corresponds a vector r. This relation is 
not uniquely reversible. Indeed, if 


z i “ ^i e 


and 


z 2 - P 2 e 


Then, from Equation 12: 


x - Pj 2 " P 2 2 , 
y = 2 Pj p 2 cos (>// 2 , 

z = 2pj p 2 sin(0 2 -0 X ) , 
r = Pi + P 2 2 * 

0 2 “ v//j - 20 . (13) 


t - AXIS 



REAL AXIS 


Figure 2— Representation of spinor Z 12 
In complex plane. 


Therefore (see Figure 2), if we symbolize the spinor s 12 by 
the arrow Zj-Zj in the complex plane or by the triangle 
Ozj z 2 ,the vector r is completely determined by the two posi- 
tive numbers p 1? p 2 and the angle 20 between the directions 
Oz j and Oz 2 . That is, the vector r is completely determined 
by the geometric form of triangle Ozj z 2 and its sense of 
circulation, but is independent of the triangle* s position in 
the plane, as given by angle 0 r (1/2) (0 t + 0 2 ). Therefore, 

0 is indeterminate and may be considered as an arbitrary func- 
tion of time. Let S* 2 be one individual spinor selected from 
the multiplicity of spinors related to the vector r. This 
multiplicity may be expressed by S 12 = S,* 2 e^, where 0 may 
denote any angle between 0 and 2tt. 


We observe that the equations of transformation, Equations 12, go over to the Levi-Civita 
transformation of Equations 7, in the case of plane coordinates. Then z = 0 and 0 2 = 0jOr 0 X + v. 
Thus, the two complex numbers are situated on a straight line through the zero-point of the complex 
plane, either on the same side or on opposite sides of the origin. The position angle of this line is 
0 which, being arbitrary, may be set at zero. Then z l and z 2 lie on the real axis and are real 


12 



numbers f and 77 , such that 


y - 2^77 . 

It is instructive to write down in spinor form some functions of the Cartesian coordinates and 
their derivatives, which are important in the theory of celestial motions. Such functions, in the 
case of undisturbed orbits, are: 

the coordinates 


xy - yx , yz - zy , zx - xz 

of the constant vector r x r of the integral of areas, 
the square of velocity, 


+ y 2 


+ z 


2 _ _ _ 1 _ 
r a » 


which is an essential part of the integral of energy. Using Equations 12, the derivatives with re- 
spect to any independent variable s (denoting derivatives by primes), and the abbreviation 


P - Zj Zj' - Zj Zj' + Z 2 z 2* “ Z 2 z 2 , 


(14) 


gives 


r ' - 2(zj Zj' + Z 2 z 2 ) + P D 2 ( z i z i /+z 2 z 2) ~ p > 


x /2 +y' 2 +z' 2 - 4r(zj' Zj' + z 2 ' z 2 ') + P 2 , 


<< 

N 

1 

N 

<< 

= 

2i 

[( Z l Z 1 

- Zj Zj') r + z 3 

, 2 .P]- 

zx' - xz' + i(xy' 

- yx' 

’) 

= 2i 

[(*!*»' ~ Z 2 Z l) 

r + z j "z 2 P J , 

zx ' - xz' - i(xy' 

- yx' 

) 

= 2i | 

_( Z 2 Z 1 _ Z J Z 2 ) 

r + Zj z 2 pj . 


(15) 


The last two of Equations 15 are each other's complex conjugates; replacing z lf z 2 , z 19 z 2 , i, and 
P by their conjugates z,, ~z 2 , z 19 z 2 , -i and P(= -P) in either equation changes it into the other. 
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Since the vector r in the xyz-space corresponds to a multiplicity of spinors S 12 , this relation 
may be fixed by imposing a simplifying condition to be fulfilled by the elements of the spinor, 
namely thatP = 0, so that Equations 15 become: 


r' - 2(zj Zj' + z 2 z 2 ') 

= 2(zj zi +z 2 z 2 ') , 


x' 2 + y' 2 + z' 2 = 

4r (zj Zj + z 2 z 2 ) , 


yz f - zy ' = 2ir 

(z, Zj' -z, z-) , 


zx' - xz' + i(xy' -yx') 

= 2ir( Z j z 2 ' -z 2 z,') , 


zx' - xz' - i(xy' -yx') 

= 2ir(z 2 z/ -z, z 2 ') . 

(16) 


But first we must prove that this does not affect the universality of vector r, i.e., that the relation 
P = 0 is independent of Equations 12 which define r as a function of spinor S 12 . 

Since 


and 


Z x ~~ p.e 


i 




z 2 “ P 2 e 


then 


P - 2i(p 1 2 </>,' + p 2 2 . 


If P = 0 and P 2 -p 1 - 2<t>, then 


Pi 'Pi + P? P 2 ~ 0 , 

Pi ~ Pi = -20' - (17) 


From Equations (17): 


Pi 




4>‘ . 


Pi 


2pi 


P' 
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If 0 = ( 1/2) (0j + 0 2 ) , an angle by which the position of the spinor in the complex plane is fixed, 
then 


0' 


Pi - P 2 


0 ' 


X 

T*' > 


with arbitrary initial values 0) — for instance 0(0) = 0. Therefore, as 0 is a well-determined 
angle connected with the movement r(s), the position angle 0 is fixed for every value of the in- 
dependent variable s. The first of Equations 17 involves the statement that the movement of the 
spinor that corresponds to the condition P = 0 may be interpreted as a motion of a system of two 
equal masses in the p, 0 -plane around the origin 0, the moment of impulses being constant and 
equal to zero. 

The next step is to write in spinor form the differential equations for Kepler T s Laws. We 
define the independent variable s by the differential relation 


dt - rds . 


Transforming the differential equations 


X + 


X 



0, 


directly into spinor equations is elementary but complicated. The transformation is simpler to 
realize starting from the known integrals of the two-body problem. These integrals give four 
algebraic relations which, by Equations 16, are: 


yz-zy - — (yz'-zy') - 2i (z y z j' - z 1 z j') - const. 


zx - xz + i(xy - yx) - y [zx ~ xz ' + i( xy ' - yx ' )] - 2i(zj z 2 ' - z 2 Zj') - const . 


1 


ZX - xz 


z - i(xy - yx) - y [zx' - xz ' - i(xy ' - yx ' )] - 2i (z 2 -z, z 2 'j - const. 


x 2 + y 2 + z 2 - y - 

x' 2 + y ' 2 + z' 2 1 


x' 2 + y ' 2 + z' 2 _2 

_ 2 r 


4r 


2 Z 1 Z 1 + Z 2 Z 2 2 


const . 


r 

4a 


(18) 
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Differentiating the first three relations of Equations 18 gives 


z i z i Zj Zj 0, 


Z 1 Z 2 “ Z 2 Z 1 “ 0 i 


0 . 


These equations let z 2 ", z^', z 2 % be expressed in terms of z{: 


z 


// 

1 





(19) 


The last of these equations can be written 


Zj z 2 z 2 Zj 0 i 


delivering the spinor form of the integral of areas 


Z 1 Z 2 ~ z 2 Z l ~ const . 

Differentiating the fourth of Equations 18 gives 


+ z 0 


-1 + Z 2 Z 2 


r 

4a 


Using Equations 19 to eliminate z 2 ", z 2 " and multiplying by z x gives 


( z i 


+ z 0 


+ z f 


r 

4a 


But ? as the bracket equals r 


z 


1 



0 , 


and likewise (multiplying by z 2 /z x ) 


z 


2 



0 . 


(20) 
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These are two linear second-order differential equations for the spinor S 12 
the spinor by its elements in vector form, 


(z j , z j j . Defining 



gives four differential equations: 

u." + CO 2 u. = 0 , (i - 1, 2, 3, 4) 

2 - J_ 

" “ 4a • 


(21) 


which may be solved like the equations for <f, 77 in the plane problem. 

Solving these equations gives the corresponding solution for the Cartesian coordinates, x, y , 
z by Equations 12. However, the same problem expressed in space coordinates is perfectly de- 
termined by six independent constants, while Equations 21 require eight initial conditions, ( 0 ), 
u.‘ ( 0 ). Two of the integral constants seem superfluous. But because of the relation P = 0 be- 
tween u. and u' , only one of the eight initial conditions may be arbitrarily chosen between certain 
limits without influencing the character of the solution. Figure 2 makes this clearer. Triangle 
Ozj z 2 , which represents the spinor, corresponds to a vector (x, y, z) only by its geometrical form 
and is independent of its position in the plane. The condition P = 0 is identical with a first-order 
differential equation determining the rotational or oscillatory motion of the triangle around the zero- 
point of the complex plane, but this motion may begin with an arbitrary initial position. 

Another question remains unanswered. In the equations 

+ CO 2 £ = 0 , T )' 1 + oo 2 Tj - 0 , 

(which may stand now for the spinor equations of an undisturbed plane planetary motion), the in- 
tegral constant of energy appears already in the equations themselves in the form of the constant 
co 2 = l/a. This is confusing, indeed. However, these equations may be considered independent, as 
each of them is a differential equation in one variable (<f or 77 ) only, with co 2 as a parameter, which 
may be different for and 77 , and is not a constant of integration. On the contrary, co 2 represents 
the magnitude of the acting force causing an oscillatory movement of £ or 77. On the assumption that 
there are two different forces co 2 for £ and co 2 for 77 , the motion of the system will be composed 
of two different oscillations with different periods for £ and 77, and this process cannot be considered 
as the Levi-Civita transformation of a motion in a Keplerian ellipse. This is only possible if the 
two parameters are identical and we consider the two equations as an inseparable system. Then 
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co 2 = k 2 (mj +m 2 )^4a is likewise proportional to the acting gravitational force. The remarkable fact 
that the integration constant, a , appears connected with this parameter, while it does not appear in 
the differential equations for the Cartesian coordinates, is due to the transformation of time t into 
pseudo-time s by ds = dt/r , as a is the mean value of r. 

THE DIFFERENTIAL EQUATIONS OF THE SPINOR OF A PERTURBED KEPLER MOTION 

The preceding sections have shown that for undisturbed planetary motion the introduction of 
spinors, together with a pseudo-time as independent variable, is advantageous. The spinor s(u i ) , 
i = 1, 2, 3, 4, where the elements u. satisfy the linear differential equations 

u." + co 2 u i - 0 , ( co 2 - constant and real) 

readily gives the Cartesian-coordinate solutions of the two-body problem from given u.. The 
solutions u. (s) are symmetric, simple, regular, and of the same form for all types of orbits. If 
the motion is in the xy-plane, we can do without spinors by using the Levi-Civita coordinates 
u 2 , where 


x - Uj 2 ~ u 2 2 , y - 2u 1 u 2 . 

Consider, now, the more complicated problem of planetary motion perturbed by one or more 
celestial bodies. Consideration of a single perturber is sufficient, as the transition to more than one 
is simple. 

Let r(t) be the heliocentric position vector of a planet. Then 

7 = --^ + p (r = | r|) (22) 


is the differential equation of its motion. The perturbation vector p(t)may be decomposed into 
three components 


p - Fr + Gr + h(? X r). (23) 

The components have the directions of the position vector r, the velocity vector?, and a normal 
vector (r x r ) perpendicular to the plane determined by r and r. The coefficients F, G, H are scalar 
functions of time t and of the heliocentric coordinates of the perturbed and perturbing bodies. 

For the purpose of regularization we could define, as before, pseudo-time s by 

ds = ydt . (24) 
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But now we follow a proposal of P. Rustaanheimo (Reference 3), who more generally defines s by 
the relation 


ds 


A -JoKdt 
T e 


dt , 



(25) 


with A, B as arbitrary constants and K(t) as an arbitrary function oft. IfA = l,B = K=0, Equa- 
tion 25 changes into Equation 24. 

Denoting derivatives with respect to s by primes, we have 


t ' 


- JL = _L Jo tKdt 


t " = 


_ J_ Jo Kdt 


[r' + rKt ' ] 


t ' 


+ Kt ' 2 . 


(26) 


Let 




Eliminating t " by Equation 26 and r by Equations 22 and 23 gives 


’■ - ( r - ?) 


7t' 2 + 




— + (K + G) t ' r' + Ht' (7x?), 


„ r ' t ' , , 

t " - — — + Kt ' 2 . 


(27) 


This is an 8th-order system composed of one vectorial and one scalar differential equation, both 
second-order. This system may be reduced to 6th-order by assigning special values to constants 
A and B. In the following, A = 1 and B = 0. 

Now let the elements of a spinor S = ( z lt z 2 ) replace the position coordinates x, y , z of the 
perturbed body. The two complex numbers z. ( j = 1, 2) will be solutions of two differential equa- 
tions, which are to be formulated. Spinor algebra permits the direct transformation of Equations 
27 into spinor form (Reference 3). The spinor equation is shown to be identical with the first of 
Equations 27; this proves its validity. 
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Assume that the elements z 1 and z 2 of a spinor are solutions of the differential equations 


z " 

j 



1 

2r 


(r' +P) 


z . 

j 


(j = 1, 2) . 


( 28 ) 


Here 


p - z y z \ ” z i z i + Z 2 Z 2 - Z 2 Z 2 


is the function previously introduced in Equation 14; f and g are two arbitrary complex functions 
of s. As P = -P, the equation for the complex conjugate z. is 


( ft ' 2 + z j + (et' -t) z ' - W( T ‘ ~ p ) ■ 


(j = l. 2) . 


(29) 


Differentiating P gives 


P ' - Zj Zj ; - z 1 Zj" + z 


2 Z 2 Z 2 Z 2 


Substituting the right-hand sides of Equations 28 and 29 for the second derivatives of the z and z 
gives 


( ft 2 + 4^) ( Zl Z > + + ( gt ' + t) j*i' z i 


r ' + p _ _ 

+ z 2 z 2 2r ( z i z i + z 2 z j) 


“ ( ft ' 2 + ^)( Z t z l + z 2 z z) - ~t) 2 


r 1 - P 
2r 


( Z 1 Z 1 + Z 2 Z 2 )] 


From Equations 12 and 15, 


r = z lZl + z 2 z 2 , r' + P - 2(z/ Zj + z 2 ' z 2 ) , r # -P = 2 (z t z/ + z 2 z 2 ') . 


Therefore, 


P' = (f - f) t /2 r , 
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and if P = const. = 0, the function f must be real. In addition, let g = h + ik, g = h - ik, where 
and k are real numbers. Then the differential Equations 28 and 29 take the simpler form 


Equation 27 gives, for the x 

X " = 


z " - ft' 2 Zj + (h + ik)t' 

~n - ft /2_ z. + (h-ik)t' 

-component of the vector T, 
(f - xt' 2 + |j + (K + G) 





z 


z 


)• 

.)■ 


x' + Ht ' (yz ' - zy ' ) . 


h 


(30) 

(31) 


(32) 


From x = zjz,-z 2 z 2 : 


x" - Z," Z, + Z 2 " Z 2 ' - Z 2 " Z 2 - Z 2 " Z 2 + 2(zj' Zj' - z 2 ' z 2 ') . 


Eliminating the z" on the right side, using Equations 30 and 31, gives 


2(z,' z,' - Z 2 ' Z 2 ') + 2(z, Zj - z 2 z 2 ) ^ft' 2 -ht' 2r) 

+ ht'[lj Zj' -i 2 z 2 +Zi Zj' -z 2 z 2 ] + ikt'fzj z/ -z 2 z 2 ' -z,!/ + z 2 z 2 '] 


From Equations 12 and 16: 


X - Zj Zj - z 2 z 2 ) 


x' - Zj Zj' - z 2 Z 2 + 2 1 Z 1 ~ Z 2 Z 2 » 


t (y z ' ~ z y ' ) 


2(z, z/ - Zj z/) 


+ 


as 


Hence 


Zj Zj' - Z, Zj' + Z 2 Z 2 ' " Z 2 z 2 ' - P “ 0- 




+ z' 


2 ( Z 1 Z 1 Z 2 Z 2 ) * 
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and 


x - X 


(at- -4 


— ht' 


x' 2 + y ' 2 + z ' 2 
‘ _ 2r 2 


+ x 




, , kt' 

- (yz - zy ) — jr- 


(33) 


Comparison of Equations 32 and 33 gives 


nc , 0 r ' x' 2 +y' 2 + z' 2 

2ft 2 - — ht' * . 

r 2r 2 


F - -r t 


ht ' + — p- = -y- + (K + G) t ' 


kt # 


Ht ' 


Equations 32 and 33 can be made identical by the following substitutions: 


k = - rH 


< h = K + G , 


1 r x 2 + y 2 + z 2 

2f = f “t? + tC k + g ) + ^ 


(34) 


These conditions are independent of the position of the coordinate axes, and are also valid if we 
consider y or z instead of x. The energy-expression 


x 2 + y 2 + z 2 1 


which is constant for undisturbed orbits but slightly variable for disturbed ones, can be used to 
rewrite the third of Equations 34 as 

2f = F + ~ (K +G) + • 

r 
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K is still an arbitrary function, it may be equated with -G. Then 


2f = F + Jj, 
r z 


h = 0 , k = - rH , 


and 


ds 


1 - f G dt 

ir e 


dt . 


Then the differential equations of the spinor are 


- ft' 2 z. 


= ikt' (z.' 2F z.j 


or 


„ Et' 2 /] 

z i 2r 2 Z j V 


F - i , 
t ' 2 + -yHt ' r 


') z i 


i rHt ' z 


As F and H are quantities of the order of the perturbing mass, the right sides of these equations 
are always small. The coefficient 


Et ' 2 _ E_ - 2 /Gdt 

2r 2 = 2 6 

is not very different from e/2, as G is of the order of the perturbations. E is not very different 
from -l/2a, if 2a is the length of the major axis of the osculating ellipse. Therefore let 

V + “ 2 z i = Qj • (35) 

where oj 2 is nearly constant, and positive in the elliptic case E < 0. The Q. are expressions of the 
order of the perturbations. 

The next section discusses the above results from another point of view. Equations 35 assume 
a form in which the coefficient of z. on the left side is a true constant. Then these equations repre- 
sent oscillations modulated by small external forces and well suited for numerical integration. 

They may be also the starting point of a new theory of general perturbations, which still has to be 
fully developed. 
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VECTORIAL REPRESENTATION OF SPINORS 


The previous section discussed the differential equations for a spinor (z 1? z 2 ) corresponding 
to the position vector (x, y, z) of a celestial body moving around the sun and subjected to per- 
turbations by one of the major planets. LetF, G, H be the components of the disturbing acceleration 
their directions are those of the position vector r, the velocity vector r, and the normal vector 
(rx r), respectively. Then the two complex elements z 2 and z 2 of the corresponding spinor are 
solutions of the differential equations 


= ft 


z. + (h + ik) t' (z.‘ - V z i) (i - 1. 2) . 


Here the real quantities f , h, k are defined by 


(36) 


, E r 

2f = F + — + — (K +G) , 

h = K + G , 

k = - rH . (37) 


Also, 


E 



(38) 


is an energy expression, which is slightly variable for perturbed motion, and is constant for un- 
disturbed motion. K is an arbitrary function of time t, and can still be specified. The derivatives 
z‘ , z", t ' are taken with respect to pseudo-time s, which is related to t by the differential equation 


ds 1 - f 1 Kdt 

dt ' r e 0 


(39) 


with initial condition s = 0 for t = t 0 . 

Giving the arbitrary function K(t) different forms will specialize these equations. Two pos- 
sibilities for K(t) are of especial interest: K = -G and K = 0. As already shown, the first of these 
special assumptions leads to 

E 

2f = F + -t, h - 0 , k - - rH 

r z 
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and to the spinor equations 


z 

J 


Et ' 2 
2r 2 


z . 

J 


1 _ 

2 


(Ft ' 2 + 


iHt ' r ') z^ 


i rHt ' z . 
j 


(j = 1, 2) 


These differential equations are of the type 


z . " + co 2 z . = Q . , 

j j ^ j 

where co 2 - -Et' 2 /2r 2 is nearly constant, and the Q. are complex functions of s and have magni- 
tudes of the order of the perturbing forces. 

Consider the special case where the perturbing planet moves in the same plane as the disturbed 
body. The problem becomes a plane one as the perturbation has no component perpendicular to the 
plane of movement, and the z. are real and identical with the Levi-Civita coordinates. Therefore, 

H k - o, and 


z." - ft' 2 Zj = 0 . (j = 1. 2) . 

This system of equations has the integral of areas zj z 2 ' - z 2 z / - const., as in the case of undis- 
turbed orbits. This is a remarkable result. In its classical form, the problem of a small body 
(m = 0) moving in the plane of the orbit of a perturbing planet (Jupiter) generally has no integral 
of motion of simple algebraic form. With the conventional treatment, the useful integral of Jacobi 
is only available in the restricted three-body problem— for a circular Jupiter -orbit. Yet here 
after the introduction of spinors (i.e., Levi-Civita coordinates), such an integral exists always , 
without any reservation. Naturally, the existence of this integral is not enough to solve the problem, 
but it helps to indicate a solution. 

Consider the second special proposal for the function K(t); i.e., K = 0. Then 

ds 1 

dt = 7 • s (°) r 0 

defines the pseudo-time s in a new form. Since the work done by Sundman, Levi-Civita, Stumpff, 
and others, this form has often been accepted for a regularizing variable in celestial mechanics. 

We have 


t ' = r 


h = G 


- rH 


x 

r 


2ft ' 


r 2 F + r r 'G + E = r 2 F + r'G + E, 




ic 2 +y 2 + 


i2 ) . i = 


+ z ' : 


2r 2 


25 



and the spinor equations (j = 1,2) 


z 

J 


+ z . 

J 




i z i ( r 


2 F + 


in 


H) + rz.' (G - irH) - 


(40) 


Equation 40 represents four equations determining the four real quantities u 3 
the spinor 


S = 



u 4 that compose 


The spinor S may be replaced by the vector (u 3 ) in the 4-dimensional space R 4 . This is not as 
illustrative as its definition by a pair of complex numbers, but it is more uniform and permits 
the elimination of the imaginary i from Equation 40. 

The vectorial form of spinors has many advantages. Consider first the case of plane motion. 
The Levi-Civita transformation, 


X j = Uj 2 - U 2 , X 2 - 2Uj U 2 , (41) 

maps the two-dimensional vector (x 3? x 2 ) onto the two-dimensional (u 19 u 2 ) plane. Differentiating 
gives 


dx 3 - 2 (u x du x - u 2 du , 
dx 2 = 2 (u 2 du 3 + u 2 du 2 ] , 


or, in vectorial form, 



( 42 ) 


with the transformation matrix 



M is orthogonal and has the norm u x 2 +u 2 2 . Integrating Equation 42 gives 



From vector and matrix algebra, 


= yz^i 2 

- 2u. 2 , 

= Ix.' ! 

1 

= 4 2 u. 2 £u.' 2 = 4r£u.' 2 


Vectors in space present difficulties in mapping R 3 onto the four -dimensional R 4 . According to 
Hurwitz, transformations of the form in Equation 42 can be made only for spaces of dimensionality 
1, 2, 4 and 8. Therefore, 



where 



is orthogonal and has the norm 2uj 2 . This leads to 


(dx ; ) = 2M(d Ui ) , 

*7 = r = Su i 2 - (43) 

2x.' 2 = 4r 2u.' 2 
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These relations correspond exactly to those of two dimensions. Explicitly: 


X 2 = 2 ( U 1 U 2 _U 3 U 4) ' 

X 3 = 2 ( U 1 U 3 +U 2 U 4) ’ 

X 4 = 0 


and 

dx l - 2 (u x duj - u 2 du 2 - u 3 du 3 + u 4 du 4 ) , 

dx 2 = 2 (u 2 duj + du 2 - u 4 du 3 - u 3 du 4 ) , 

dx 3 = 2 (u 3 du x + u 4 du 2 + Uj du 3 + u 2 du 4 ) , 

dx 4 = 2(u 4 duj - u 3 du 2 + u 2 du 3 - Uj du 4 ) = 0- 

The expression on the right side of the last equation is not a complete differential. Equating it 
to zero gives the same relation, namely P = 0— previously assumed as an accessory condition 
to make the relation between vector and spinor unique. Here, this condition emerges naturally. 

We can map the three-dimensional x.- space onto the four-dimensional u. -space in the following 
manner: from the four-dimensional x.-space take out the three-dimensional section x 4 = 0; then 
the corresponding u. -space will be specialized by the differential condition P' = 0. 

We can transform the differential equations in x . directly into those in the corresponding u. - 
coordinates— for instance, by transforming the forces Pj (xj • • ■ x 3 ) into corresponding forces 
S i (uj * * * u 4 ) , where the work done by the forces must be independent of the coordinate system, i.e. 


3 4 

p i dx i = s i du i • 

i= 1 i = 1 

A most interesting study based on the above mentod is presented by Kustaanheimo and Stiefel 
(Reference 2). 

Considering Equations 43, write Equation 40 in the form 


z 

J 



1 

2 


z. (Fr 2 + irr ' 


H) 


+ z / r(G- irH) 


rR. . 

j 
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Let 


Ri - Sj + iS 2 , R 2 - S 3 + iS 4 


Then the spinor equation can be rewritten as four equations involving only real quantities 


U I + 7'(t~ lu l 2 ) u i “ rS i 


(i = 1, 2, 3, 4) 


( 44 ) 


where 


= Xu 2 , 

j 


2 2 u . u * 
j j 


S l(u,, U 2 ) 
S 2 ( U J, U 2 ) 


2 Uj Fr ~ 2 u 2 r'H + u i^ + u 2 


~2 u 2 Fr + ~2 u j r 1 H + u 2 G - u' rH 


s 3 - s,(u 3 , u 4 ) , 

S 4 = S 2 ( U 3- u 4)' 

These equations have the form of Equation 35: 


U . + CO U . 


s. 


where go 2 is nearly constant. Instead of s , introduce another independent variable q such that 
Equations 44 take the form 


d 2 u . 

J 2“ + 0 J 2 U = Q. , (i " 1, 2, 3, 4) 

dq 1 

with constant J 2 . This is the case of simple harmonic oscillation modulated by small external 
forces Q. . Indicating derivatives with respect to q by bars gives 


du 
dq ’ 


s 


ds 

dq 


Then 


, _ du _u_ 

u " dT ' s • 


d 2 u _ u us 
ds 2 ■ -s 2 " s 3 
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Then, instead of Equation 44. 


If 




r s 2 S. + u . ■==" 


v (4 -*=,’) 


is a constant, then, 


s 2 - 2 | TO) 2 + S u . 2 j , 


(45) 


(46) 


and Equation 45 gives 


u. + or U. 


= 2 rSj (ru> 2 + Z.U. 2 } 


__ s 


From Equation 46, with constant co 2 , 


S S - TO) 


+ 2 2u • u . 


2 Su. . + w 2 u 


Introducing Equation 45 gives, finally, 


1 _ = 

2 s s 


r s 2 !S. u . 

) 3 



(47) 


This may be written 


t(t s 2 - lu f) - r s 2 £S. Uj 
which, with the aid of Equation 46, becomes 


.2 - 


r s 2 2S. u. 
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Therefore, instead of Equation 45: 


u.+oj 2 u. - r S. s 2 + u . — - IS. u. 

J J 


= 2 (r co 2 + Iu/)(r Si ^XS.U.) - Qi 


(48) 


This is a system of equations of the required form. To regain the time t from the pseudo-time q, 
put 


dt _ dt ds 

dq ~ ds dq 


r /2(r<u 2 + IWf) 


The square root is not very different from unity. For undisturbed motion, co 2 = l/4a, and 


X u . 

j 


2 



Vi 

4r 



with 


V 2 



From this relation: 


X u 




i.e. 


Xu 2 

j 


2 “ 4a 


With tco 2 = x / 4a, it follows that, for undisturbed motion, 


dq ” 1^2 (tco 2 + 2u. 2 ) - 1 


As the constant oj 2 can be arbitrarily chosen, it will always be suitable to put co 2 = 1/4 a, if a is 
the semi-major axis of an osculating or intermediary Keplerian ellipse. 
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This method of using spinors in their vectorial form to solve the problem of computing disturbed 
orbits has been successfully applied to minor planets at the Institute of Applied Mathematics at 
Zurich, Switzerland (Reference 4). Difficulties arise in the case of orbits with parabolic or near- 
parabolic character. Then oP - 0 and on the right side of Equations 48 co 2 is in the denominator. 

This drawback disappears only if IS.u. goes to zero with oj 2 . The vanishing of Is u. occurs for 
all cases of undisturbed motion, as s = 1, and therefore ¥ = 0. But in general, if co = 0, Equation 
46 shows that s = 1 only when Su. 2 = const. = 1/2. This agrees with what was said about Equation 
8. Thus, with spinor coordinates and with pseudo-time q instead of t, the expression for the gen- 
eralized kinetic energy is constant in the parabolic case; therefore, the work done by the acting 
forces (defined as the variation of the kinetic energy with time) 

£s i*i = 0 • 

The above drawback is not critical; there are many ways to overcome it. For example, can 

be developed in a series 


SS.u. - a^ap + a 2 aj4 + a 2 Ctj6 + ■•• , 

that converges rapidly for small oP and has w 2 as a factor. 

The next section discusses some ideas about regularizing dynamical processes in general. 
The difficulties with highly eccentric orbits are found avoidable from the beginning. 


GENERAL CONSIDERATIONS ON REGULARIZATION IN DYNAMICAL PROBLEMS 

The preceding section presented a survey of the use of spinors in celestial mechanics. Intro- 
ducing these symbols, together with a certain pseudo-time as independent variable, made the equa- 
tions of celestial movements become simple, linear, and with regular solutions even in the case of 
collisions; which would not occur with Cartesian space coordinates and time as variables. Before 
going more deeply into the questions of linearization and regularization, let us review the useful 
attributes of spinors. 

Spinors are most obviously useful for analyzing undisturbed motions of a body around its 
central mass point. The equations of motion of the body with respect to the central mass are 

. x .. y z 

x + — r 0 , y + — t- 0 , z + — r- - 0 

r r 6 i " 3 

if the units of time and mass are so chosen that the constant expression k 2 (m x + m 2 ) is equal to 
unity and therefore disappears from the equations. The distance r satisfies the second-order 
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differential equation 


where a 2 is a constant equal to 1/a, a being the major semi-axis of the conic section orbit de- 
scribed by the body. These equations can be transformed into a linear form if pseudo-time s in- 
stead of time t is used as independent variable, such that 


dt - r ds . (s - 0 for t - 1 0 ) 

Then the differential equations for r and x take the form 


r" +a 2 r - 1 


x" + a 2 x - T 


1 

7 (x r 


xr") - const . 


The equations for q r r, x, y, z respectively, have the common form 


q'" + a 2 q' - 0 


1 - r„ 


^ ds 


( 49 ) 


The solutions of these linear differential equations may be written 


q(s) - q 0 + CjSqg' + C 2 sq 0 " 


( 50 ) 


Here q 0 , q 0 ', q 0 " are the values of q, q', q" at the initial time t = t 0 , and 




^ 2 ) • 


belong to the series of c -functions 


c Q ° cos X 


»in X 


1 ~ cos X 


X - sin X 


( 51 ) 


which are derived from the generating function c 0 = cos \ by 




cos X (ds ) n 
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In particular, 


r ( s ) - r 0 + c lS r 0 ' + c 2 s 2 r Q " , 

and the connection between intermediate time r = t - 1 0 and pseudo-time s is given by the "main 
equation" 


T 


f 

•'O 


r( s )ds 



+ c 2 s ^ r o") ds ; 


x.e. 


r - sr 0 + c 2 s 2 r 0 ' + c 3 s 3 r 0 '' , 

if we use the integral relations for the c -functions, 


f 

Jo 


ds 


~ n + 1 


> n + 1 


Equation 50 (the solution of Equation 49) is regular with respect to pseudo-time s . Indeed, in the 
case of rectilinear motion the velocity V(s) = fx ' 2 +y' 2 + z ' 2 is zero at the time of impact, while 
V (t ) = /x 2 + y 2 + z 2 ” approaches infinity. On the other hand, there remains a certain singularity of 
position. For when the falling body reaches the central mass point, the pseudo-velocity goes to 
zero and the direction of motion reverses. The space trajectory therefore has a cusp at this 
moment. In any rotating system the trajectory will appear as an ejection orbit. This singularity 
is overcome by using the elements u. of a spinor instead of the Euclidean coordinates, for then: 

1. The differential equations take the form of Equation 21: 

u + o) 2 u. ~ 0 , co 2 - 4^" - a 2 , 


and are satisfied by general solutions analogous to Equations 10: 


u . 


+ c . su . \ 


c c (a; 2 s 2 ) 

n n \ ' 


which are noticeably simpler than those of the ordinary rectangular coordinates. 
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2. The singularity of position in the case of impact disappears. In terms of the Levi-Civita 
coordinates 77, the body passes the point of impact in a straight line with finite constant velocity 
/ 0 . 


This theory, as proposed by Kustaanheimo and Stiefel, was then extended to the more compli- 
cated case of perturbed orbits. The perturbed differential equations of the spinor elements u. 
take the form 


u." + cj 2 u. - Q. (s) , (52) 

the (s) being small quantities of the order of the perturbing masses. The "frequencies” o? are 
no longer constant, but are slightly variable functions of s. To avoid this inconvenience, we trans- 
formed Equation 52 once more so that oo 2 became constant. This was made possible by introducing 
instead of s a new pseudo-time q. If derivatives with respect to q are indicated by bars, 


du 

dq 


the differential equations take the form 


u i + " 2u i - S.(q) , 

the s. ( q ) being small quantities of the order of the perturbations, and the frequency & now being 
constant. Parabolic orbits present some difficulties. The s. contain a term of the form cr/cJ 2 
which, in the parabolic case, co 2 = 0, seems to have the indeterminate value 0/0, but approaches 
a finite-limit value as co approaches zero; it may be developed for small values of co 2 into a rapidly 
converging power series. But this would necessitate the use of different computing methods for 
different types of orbits. This drawback is due to the introduction of q instead of s; fortunately, 
there is another way of making the coefficient of u. constant, without changing the variables again 
(see Reference 6, page 573 ff.). 

Consider the perturbed differential equation of r(s), 

r'" + a 2 r‘ = R( s ) = rf(s), a 2 “ ~ , 

where f (s) is a well-known function of the coordinates of the bodies that participate in the motion. 
Then 


da 2 

ds 


- J (r'" + a 2 rj - - f (s) 
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is not zero as it would be in the undisturbed ease, f = 0. Integrating this equation gives 


- - { 


f ( s ) ds 


where a 0 2 is a real constant of integration. Therefore, 


r"' + a 0 2 r' 


rf(s) + r 


fds ■ ds f r J fds l - h '( s ) 


or, after integration, 


r" + a 0 2 r - h(s) - 1 + k(s) , 

where k(s) is of the order of the perturbing masses and can be developed into a power series or a 
trigonometric -function series. Using the c -functions permits the special solutions of this equation 
to take the form 


r(s) - r 0 + CjS r 0 ' + c 2 s 2 r 0 " + (c 3 s 3 k„' + c 4 s 4 k Q " + • • • ) , 

the terms in brackets being small, and of the order of the perturbations. Similar solutions can be 
found for the rectangular coordinates x, y, z as functions of s , and likewise for the oblique-angled 
coordinates F, G, H defined by the relation 

r = Fr Q + Gr 0 + H(r 0 x r 0 ) , (53) 

where the directions are those of r 0 , the position vector at t = t 0 , the velocity vector 7 0 , and the 
vector perpendicular to the plane of r 0 and r 0 , respectively. (Functions F, g, h in Equation 53 
are different from those in Equation 23, which are the components of the acceleration p of the 
perturbing vector.) The intermediate time r = t-t 0 can be calculated as a function of s from 

si- 0 + c 2 s 2 r 0 ’ + c 3 s 3 r 0 " + (c„s 4 kr fl ' + c 5 s s k 0 '' + • • •) , 

which is the main equation of the problem. 

This elegant solution of the special problem of perturbed planetary motions can easily be ap- 
plied to the spinor equations of perturbed motion. In fact, the method can be extended to any 
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general dynamical problem whose differential equations can be transformed into linear differential 
equations by introducing new dependent and independent variables; or, the problem may have this 
form from the beginning (see Reference 7). Assume that a dynamical process p(s) may be described 
by the equation 


+ a 2 p( n ~ 2 ) + **• + a n _ 1 p' + a n p - 0 , (54) 

the being constants and s being either time or any pseudo-time. Then, if p x , p 2 , * * • , p n are n 
independent, simple, particular solutions of the problem, a generating linear function of the p. 
may be written: 


4>(s) ~ a l p 1 + a 2 p 2 + ■•• + a n p n , 

the coefficients a. being constant or zero. From this generating function, a system of auxiliary 
functions ^(s), <£ 2 (s), ••• 4> v (s), ••• may be derived by putting 

<£ v ( s ) = ~T [ *** f 0(s)(ds) v , v - 1, 2, ••• (55) 

S Jo Jo 

As 4> is also a particular solution of Equation 54, it satisfies the equation 

0( n ) + a ^(n-l) + a ^0(n“2) + ••• + a n _ 1 = 0 . (56) 

Any special solution of this equation may be given by n initial values for t t o? 

■ ^o"’ ■ ■ ■ • 4 > o n ~ 1) ■ 


Determine these values so that 


4>o = 1 


‘V + a A “ 0 • 

4>o + a i ( t>o + a 2^0 = 0 


*0 


(n-1) 


a A' 


(n-2) + 


+ 3 n-l ^0 


(57) 
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Then, after consecutive integration of Equation 56, 


<£(nrl) + Si 0<n-2) + ... + a R _ 1 0 + 3 n f 

Jo 


4> ds - 0 , 


^j(n-2) 4 a i (^( n “'3) + + a n -l ds + a n ds 2 - 0 


i: 


.f-r 


4> 1 + a x 4> + a 2 I <£ds + + a 


^ds"’ 1 = 0 


“ J> * -II 


(^ + a. I <[> ds + a 2 I I ^ds 2 + *** + a 


.f-f 


<£ds n - 1 


(58) 


The last of Equations 58 can be written (if the integrals are eliminated by Equation 55): 


<t> + BjS^j + a 2 s 2 <£ 2 + ••• + a n s n 4> n ~ 1 


Integrating v times between 0 and s, and considering the initial conditions, Equations 57, gives 


S^ v + a 1 s 1/+1 <£ v + 1 + a 2 s v + 2 4> v + 2 + 


• a n^ V+n = — 


or, dividing by s v , 


a I s ‘fiv+l + a 2 s2< ^v + 2 + ••• + a n s n t^ v+n 


(59) 


Substituting the left side of Equation 59 for each 1/v! -type factor in the Taylor series: 


p(s) - P 0 + i'j P 0 ' S + Ji p 0 "s 2 + 31 p 0 "'s 3 
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reduces tins infinite power series to the finite form. 


p(s) - Po + s^ 1 p 0 ' + s 2 ^> 2 (p 0 " + a t p 0 ') +••• 

+ sE1< ^n(pO (n) + a i PO^'^ + •** + a n-lPo') 2 

because higher-order terms on the ri^it side are zero. 

This method can be shown to apply to many dynamical problems, for instance, the three fol- 
lowing examples. 

1. Equation 

p'" + a 2 p' = 0 , ( a i = °> a 2 = a Va 3 = °) 

of the two-body problem, the generating function being <f> = cos as, which fulfills the conditions 

~ 1 » <V * a 1^0 - 0 , 4 >” +• a r <£ 0 ' + a 2 g3 0 = 0 

and the recurrence formula 

d> t <x 2 s 2 d> — — p . 

The <p v are identical with the c -functions. 

2. The differential equation of a periodic oscillation with two frequencies, 

p a t cos ^a x s + AJ +■ a 2 cos (a 2 s + A 2 j . 


This function is the solution of the linear and homogeneous differential equation 




which obeys the conditions 


V + a 1^0 :::: 0 • 


* u 


- (VW) . 
0 , 


f ^ i»i .-A ^ ^ i-i ■/■, ' ** 4- ■ ■£, ‘ /■■} "" fl 

r 0 a 1 r 0 a 2 v 0 a 3 r 0 a 4 U ’ 


* Ks 
+ 0 


a 4 + a.“ a-“ * a 4 


3. The differential equation of a damped vibration, 


+ 2kp /? + (k“ + a 2 ) 


In this case the generating function is 


|cos a s " ” sin a s j 


"k s f 


from which the auxiliary functions 


,1: 

'r 



■±<ds r 


C£in be derived. Then 


h + 2k s 4 \, + j + (k 2 + a 2 ) s 2 ^ 


eliminates the reciprocal factorials from the Taylor development of p(s) and produces a closed 
formula. In these last two examples s is the time t - t 0 itself, but we can imagine cases in which 
the simple form p(s) appears only after a convenient transformation of the coordinates. 

The above considerations help to solve the two-body problem and may apply to more compli- 
cated problems. They have not been applied so far, but it would be very well worth while if they 
were. This method may be qualified for finding, for instance, closed mathematical expressions 
for certain periodic orbits in the restricted three-body problem. 

The author of this paper has applied a similar method to solve Hill’s problem of the periodic 
inequality of the motion of the moon or M variation" (Reference 6, pages 578-588). There are other 
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ways of using these ideas. For instance, the calculation of intermediate orbits of the plantes or 
the moon; this includes calculating the osculating Keplerian orbit and some of the more important 
periodic terms of the perturbations. And doubtless there are many more. 
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